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We report an experimental implementation of the minimal Deutsch algorithm in an optical setting. 
In this version, a redundancy is removed from the most famous form of the algorithm. The original 
version involves manipulation of two qubits, while in its minimal version, only one qubit is used. 
Our qubit is encoded in the transversal spatial modes of a spontaneous parametric down-converted 
signal photon, with the aid of a double slit, with the idler photon playing a crucial role in creating 
a heralded single photon source. A spatial light modulator (SLM) is programmed to physically 
generate one-bit functions necessary to implement the algorithm's minimal version, which shows 
that the SLM can be used in future implementations of quantum protocols. 



I. INTRODUCTION 

Quantum computation has emerged in the past 
decades as a potentially powerful tool to solve problems 
more efficiently than its classical counterpart. One sim- 
ple example is determining whether a coin is fair (heads 
on one side, tails on the other) or fake (heads or tails 
on both sides). This is one version of the Deutsch prob- 
lem [![, who himself showed that, when exploring quan- 
tum state superposition, only one examination step is 
necessary, while classically the solution requires individ- 
ual examination of both sides. Although the algorithm 
originally proposed by Deutsch involves manipulation of 
two-qubit states, there is also a "minimal" version of it, 
in the sense that just one qubit is manipulated @. In 
this sense, minimal Deutsch algorithm can be considered 
the most basic and simple quantum computation. 

Many physical systems have shown to be useful for im- 
plementing quantum computation, such as nuclear mag- 
netic resonance 0, 0| , trapped ions [f| , optical cavities 
[(| , Joscphson junctions Q , and photons @, . In par- 
ticular, the original Deutsch algorithm and its generaliza- 
tion, the Deutsch- Jozsa algorithm fiol ] , were implemented 
using photons [1, GH], and its minimal version was also 
implemented on quantum dots [l2j . 

Spontaneous parametric down conversion (SPDC) (l3j 
is a natural source of correlated photon pairs, with the 
additional advantage of having many degrees of freedom 
that can be considered quantum systems, such as polar- 
ization [H,[llj], transversal [lU, longitudinal HE Ell; and 
orbital angular [lU, [l9| momenta. On the other hand, a 
spatial light modulator (SLM) can be used to perform 
state control [201 . It has been used for tomographing 
polarization [2ll] and transverse momenta (22j states, for 
measuring Bell inequality violations in orbital momenta 
of SPDC photon pairs [23( , and optical quantum algo- 
rithm simulation |24| . 

In this work we report an experimental implementa- 



tion of the minimal Deutsch algorithm in an optical set- 
ting. We use a double slit to encode logical qubits (in 
the sense of the {|0) , |1)} logical base of a qubit space) 
in the transversal spatial modes of photons generated in 
a SPDC process [15|. A SLM is employed as a fundamen- 
tal part of our experimental setup. It has the function of 
simulating the "fair coin" or the "fake coin" in the opti- 
cal setup. Before really applying the quantum algorithm, 
we need a calibration process, which can also be under- 
stood as a proof of principle of the algorithm, since one 
uses many caries of the oracle. After such calibration, 
the apparatus is ready for running the real Deutsch algo- 
rithm: With only one "examination step" , answer which 
type of "coin" we have, with probability larger than 1/2. 
The paper is organized as follows: In Sec. II we review 
the theoretical description of the minimal Deutsch algo- 
rithm. The experimental setup is presented in Sec. III. 
The results are shown in Sec. IV. Discussion is made in 
Sec. V and conclusions arc outlined in Sec. VI. 



II. MINIMAL DEUTSCH ALGORITHM 

Quantum parallelism allows quantum systems to eval- 
uate a function f(x) for many different values of x si- 
multaneously. The Deutsch algorithm is a good example 
of how to explore quantum parallelism to answer a clas- 
sical question; explicitly, to solve the Deutsch problem 
evaluating the function only once. 

Consider an oracle that can answer one-bit questions 
with one bit answers described by a deterministic func- 
tion / : {0,1} — > {0,1}. This function is called bal- 
anced if /(0) ^ /(l), otherwise the function is constant 
[/(0) = /(l)]. The Deutsch problem consists in deter- 
mining whether a given unknown function is balanced or 
constant. For a classical algorithm to answer that with 
certainty it requires the oracle to be asked twice, that is, 
asking the value of / on and 1. On the other hand, 
the Deutsch algorithm requires only a single query, us- 
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ing quantum parallelism, to reduce the minimal resource 
required. 

In the minimal version , the oracle's behavior is en- 
coded in a unitary operation Uf to be applied on a well- 
chosen input state, depending on the function /. For 
the computational basis one has Uf \x) = (— l)f( x > \x) 
as output, for x = 0, 1. Now, if the superposition state 
^(|0) + |1)) is used as input, both questions are asked 
at the same time, with the output state being: 

U) = L r. l0) + Il) = |Q) + (-!)/(!) |1) 

Disregarding global phases, one has: 
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Note that the first answer is orthogonal to the second, 
so we can make a projective measurement in the ba- 
sis {|+) , |— )} and find out if the function is balanced 
(|— )) or constant (|+)). The unitary operations are im- 
plemented by the oracle. Four possible maps L7y are 
generated and we label them with the values of i = /(0) 
and j = /(l): 
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/(0) = /(l) = 0, 



/(0) = 0^/(1) = 1, 



/(0) = 1 /(l) = 0, 



/(0) =/(!) = !. 



(3) 



III. EXPERIMENTAL SETUP 

The scheme of the experimental setup is shown in 
Fig.l. A 50 mW Hc-Cd laser operating at A = 325 nm 
is used to pump a 2-mm-thick lithium iodate crystal and 
generate, by type I SPDC, degenerate non-collinear pho- 
ton pairs. Signal and idler (A s = \ = 650 nm) beams 
pass through a A/2 plate (half- wave plate), before they 
cross a double-slit placed at a distance of 250 mm from 
the crystal. The double-slit plane (xy plane) is aligned 
perpendicular to the plane defined by the pump laser and 
the down-converted beams (yz plane), with the small di- 
mension of the slits parallel to the x direction. The slits 
are 2a = 100 fim wide and have a separation of 2c? = 250 
/im. The lens L\ is used to generate photon pairs in en- 
tangled transversal path states (25[. A natural question 
is why to use down-converted biphotons for implement- 
ing the one-qubit Deutsch algorithm. In our case, one 
can consider the selective detection of the idler photon 
as part of the heralded preparation of the signal photon 
state. 
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FIG. 1. (Color online) Experimental setup scheme for mini- 
mal Deutsch algorithm implementation. The L\ lens focuses 
the pump beam in the double-slit plane; lenses L a \ and Li are 
used to detect the signal and idler beams at Fourier plane, 
while the L 3 2 lens is used to project the double slit images 
in the detector. A half-wave plate is placed right after the 
crystal and a polarizer P s is positioned in front of detector 
1. CNC denotes coincidence counter and SLM denots spatial 
light modulator. 



Our experimental setup is arranged in a way that a 
photon passing through the inferior slit of the double 
slit corresponds to state |0), while a photon that passes 
through the superior slit corresponds to state |1). The 
SLM after the double slit, together with the A/2 plate 
and the polarizer P s , works as the quantum oracle. The 
map of the oracle function is constant (Uqo or Uu) if the 
phases added by the SLM are equal, or balanced (Uqi or 
U\o) if the phase difference is tt. Once the pump beam is 
focused at the double slit plane, a Bell state cre- 
ated by the twin photons in the slit path states [2^, [26| . 
A dichroic mirror placed just after the crystal removes 
the pump beam and transmits signal and idler beams. 
The trigger photon (idler) also passes through a double 
slit and is reflected by the SLM, but without the polarizer 
at its path, no phase change or amplitude variation in its 
state occurs due to the manner in which the SLM works. 
In Fig. [5] and Fig. [3J we show, experimentally, that we 
are able to introduce spatial phase changes at the signal- 
photon path state while preserving the state amplitudes 
by using the SLM. However, there are many other maps 
that the SLM can implement in a double slit-qubit. The 
SLM maps can be described by: 



Uslm 



Aoe^ 






Aie^ 1 



(4) 



where Ak is the attenuation and <f>k is the phase applied 
in the photon state \k). The SLM maps in this kind of 
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setup are diagonal because it cannot exchange photon 
population between slits, i.e., an operation like |0(1)) i— > 
A |0) + Ai |1) cannot be done, if Ak ^ 0. Instances 
of these SLM maps were implemented in Refs. [U [22| , 
while a general one can be made through a calibration 
described by Moreno et at [27I ]. 

Single slits with 100 /xm width are placed in front of 
each detectors. Their planes (xyi^ s planes) are aligned 
perpendicular to the propagation direction of the idler 
and signal beams (zj iS direction), respectively. The small 
dimension of each slit is parallel to the corresponding x 
direction. The SLM used is a Holoeye Photonics LC-R 
2500, which has a 1024 x 768 pixel resolution (each pixel 
consists of a 19 x 19 /mi square) and it is controlled by 
a computer. Signal and idler beams are focused on the 
detectors with a microscope objective lens (not shown in 
Fig- HI)- Two interference filters, centered at 650 nm and 
10 nm FWHM bandwidth, are kept before the objective 
lenses. Pulses from the detectors are sent to a photon- 
counter and a coincidence detection setup with a 5.0 ns 
resolving time. 



IV. EXPERIMENTAL RESULTS 

First of all, we must be able to implement the maps E/y 
e {0, 1}) of Eq. © . It was shown H3] that a SLM 
plus the input and output polarizers can be properly cal- 
ibrated to obtain this goal (A = A\ and <f>k = 0,n). The 
liquid crystal display of the SLM is divided in two regions, 
each region adding a phase to the photon path state, |0) 
or |1), defined by the slits (inset of Fig. [T]). A SLM gray 
level is associated with each region of the display. Pre- 
determined gray levels, along with correct half- wave plate 
and polarizer angles, introduce a relative phase between 
the photon path states without relative amplitude atten- 
uation. The evolution maps Uij are implemented when 
the correct phase differences (0 or n) with no amplitude 
attenuation are introduced in the photon path states by 
the SLM. 

Figure [5] shows two double-slit interference patterns, 
measured in coincidence counts. Both are obtained by 
scanning the signal beam detector, with steps of 40 /xm, 
while maintaining the idler beam detector fixed in the 
position that corresponds to an interference pattern max- 
imum (xi = 0). In the closed squares pattern, we have 
used the same SLM gray level for the |0) and |1) signal 
photon path states. But in the open squares pattern we 
have inserted a relative phase between the photon states, 
through an appropriate choice of SLM gray levels for each 
one. Fitting the two interference graphs, we measure a 
relative phase of A<fi = 3.25 ± 0.03. To obtain the in- 
terference pattern in both detectors, we have used the 
Lj and L s i lenses in our experimental configuration. It 
is important to note that the phase in the idler photon 
path states is not affected by the SLM gray level, because 
there is no polarizer in front of this detector [27j . 

In Fig. [3] we have the coincidence double-slit image, 
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FIG. 2. Coincidence double-slit interference patterns. Closed 
squares show the interference pattern when the SLM gray lev- 
els are the same in both slit aperture directions that define 
the signal photon paths states. Open squares show the inter- 
ference pattern when there are different gray levels producing 
a relative phase of ir, between the signal states. Idler path 
phase, in both measurements, is not changed by the SLM. 
The idler detector is kept fixed at Xj = 0, while signal detec- 
tor is scanned in steps of 40 /jm and the detection time is 60 
seconds. Lenses Li and L s i were used. 



for the signal beam. This result is obtained, when the 
L S 2 lens is used in the signal beam and the Li is used 
in the idler beam. Here the experimental setup is such 
that the peak in the x s displacement negative region is 
associated with the inferior slit, i.e., with the |0) signal 
photonic state, while the peak in the positive displace- 
ment region is associated with the |1) photonic state. 
Once again wc kept the idler detector fixed at X{ = 0, 
i.e., its interference pattern maximum and scan the sig- 
nal beam detector, with steps of 20/im and acquisition 
time of 20s. Signal detector is placed at the image plane 
while idler detector is at the Fourier plane. The SLM 
gray levels were the same as the ones used to obtain the 
open squares interference pattern shown in Fig. f2J Wc 
can infer that the SLM gray levels used to obtain phase it 
do not attenuate the state amplitude by calculating the 
areas under the peaks corresponding to each state. In 
the curve shown, we have an area of (77 ± 2) arb. units, 
for the |0) signal photon state and (72 ± 2) arb. units, 
for the |1) state. Therefore, with the curves shown in the 
open squares of Fig. [5] and Fig. [31 wc can implement the 
Uqi map. To implement each Uj a proper gray level was 
chosen in each slit, corresponding to add a 7r phase for 
the value 1 and no extra phase for the value 0. 

To implement the minimal Deutsch algorithm we must 
create the state |+) for the signal photon and, after the 
oracle, measure it in the base {|+),|— )}■ By detect- 
ing at the Fourier plane and at the origin of the inter- 
ference pattern, we are able to implement experimen- 
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TABLE I. Experimental results for all maps possibilities. 
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FIG. 3. Signal double-slit image, measured in coincidence 
counts. The result is recorded with the idler detector fixed 
at Xj = 0, while signal detector is scanned in the x direction, 
when the same gray level of open squares pattern used in 
Fig.[5]is applied and the detection time is 20 seconds. Lenses 
Li (the same lens used in Fig. [2} and L S 2 were used. Signal 
detector is placed at the image plane while idler detector is 
at the Fourier plane. 



tally the detection projector |+) (+| [26j. Using SPDC, 
a double-slit, and by focusing the pump beam at the 
double-slit's plane, we prepare the Bell state \il> s .i) = 

l^ + > =^(l°- 1 *) + l 1 .°i» = ^(l+.+i)-|-.-i»S- If 
we detect the idler photon using the projection operator 
|+) (+| wc project the signal photon state in the state we 
need. Hence, in our setting, an oracle query corresponds 
to a detection of the idler photon at the center of the 
pattern, being the instance where the logical qubit, the 
signal photon, is prepared in the appropriate state. The 
oracle's answer is then provided by the detection or no- 
detection of the signal photon. Note, moreover, that the 
idler measurement, and hence logical qubit state prepa- 
ration, is done after one of the maps Uij is applied on the 
qubit. But since these operations commute, this does not 
affect the final statistics. 

The experimental results for all map possibilities are 
shown in Table HI The measured photon coincidences were 
obtained with signal and idler detectors kept fixed at 
x s = 0, and Xi = (center of the interference pattern), 
respectively, at the Fourier plane. Measurements were 
taken in 1000 s. The data show a clear difference of 
behavior between constant and balanced functions. 



V. DISCUSSION 

In the previous section we presented the result of the 
algorithm using many oracle queries, summarized in Ta- 
ble HI which allows, as expected, perfect discrimination 
between constant and balanced functions (high and low 
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coincidence counts, respectively). The Deutsch algo- 
rithm is, however, about the optimization of the Deutsch 
problem with regard to the number of queries to the or- 
acle; and here we discuss what our experimental setup 
tells us when only one such query is allowed. 

However, the meaning of this implementation regard- 
ing individual events is more subtle, due, among other 
reasons, to the fact that an idler photon detection at 
position x corresponds to the preparation of a state 
7j7|(|0) + e 1 ^^ |1)) for the signal photon, where (f)(x) is 
a function depending on the detector position and other 
features of the setting (see Refs. [25l. [2a| ) . Strictly speak- 
ing, the measurement is not of the von Neumman type, 
on the states {|+) , |— )}, but a positive operator valued 
measurement (POVM). This means that an infinitesimal 
detector in x — measures an operator proportional to 
|+) (+|, while a finite detector with opening d (that is, 
the width of the slit placed in front the detector) mea- 
sures a positive operator that is a weighted sum of all 
projectors from x = — | to x = |. Moreover, a signal no- 
detection event does not necessarily correspond to the 
state |— ). Nevertheless, we can use the results shown in 
Fig. [2] to predict the chances of having a detection (and 
hence, also no detection) event when the function is con- 
stant or balanced, by computing the area under the cor- 
responding (normalized) curve over an interval centered 
at the origin with the same width as the detector. 

From the graphs themselves, and also from Table HI wc 
see that the chances of detecting a photon when the func- 
tion is constant is much higher then when it is balanced, 
so indeed a single oracle query gives us some informa- 
tion about the function. But, in practice, none of the 
events tell us definitely which type of function we have. 
A detection can also be associated with a balanced func- 
tion, although it is rare, due to detector width. A no- 
detection event, on the other hand, might be related to 
one of three distinct situations: the function is balanced 
so, as we want, there is a very low probability of detect- 
ing a photon at the origin; or the function is constant, 
but the detector failed; the function is constant, but the 
photon hit the Fourier plane at another point away from 
the detector. Due to this third situation, even in a per- 
fect experimental setting a no-detection event can also 
correspond to a constant function. 

To understand better the quantum advantage within 
this implementation we can consider a scenario where one 
of these functions is given to us, with equal probability, 
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and we must bet on constant or balanced with only one 
oracle query. Classically, we cannot do better than a fifty- 
fifty guess, but with this implementation we can. Indeed, 
from the discussion in the last paragraph, we already see 
that a detection is more likely to be associated with a 
constant function, while no detection is more likely to be 
associated with a balanced one. 

To see this in a quantitative manner, let S be the event 
where the function's type is correctly guessed. We can 
write the probability for this event as: 

P{S)= J2 p (f = WP(S\f = ij), (5) 

i,j=0,l 

where P(f = ij) = 1/4 is the probability of having the 
function /(0) = i,f(l) = j while P(S\f = ij) is the 
probability of success given that function. We denote 
by Pij the probability for a photon to hit the Fourier 
plane at a point covered by the detector, given that the 
function implemented is / = ij. That is, pij is just the 
area under the normalized curve of Fig. [5] corresponding 
to the function ij, in an interval around the origin with 
the detector's width. 

Now, if the function is the constant 00, we succeed in 
our guess if we detect a photon. This will take place 
with probability r]poo, where r\ is the detector efficiency. 
Similarly, for / = 11, we have P(S\f = 11) = r\p\\. For 
the function 01, on the other hand, we succeed if there is 
no detection, which has probability (1 — poi) + (1 — V)Poi- 
The first term in the sum corresponds to the case where 
the photon goes to a point away from the detector, while 
in the second the photon hit the detector but the detector 
fails. Of course, we have also P(S\f = 10) = (1 — pio) + 
(1 — T))piq. Since the experiment is designed in a such a 
way that poo ~ Pu = Pc and Poi ~ Pio = Pb, we have 
finally: 

P{S) = \[l + r,{p c - Pb )]. (6) 

Of course, for any detector width, the probability is 
just 1/2 for r\ = 0, since we do not gain any information 
about the function. It then grows linearly with r\ and, 
for the detector width we have used (100/mi), it goes to 
a maximum of 0.55. 

We could also vary the detector's size by changing the 
width of the slit placed in front of it, to maximize the 
right function choice. Figure [4] shows the success prob- 
ability when the slit size is changed, and the detector is 
considered perfect. The best detector slit size is 260 pm 
and the corresponding success probability is 0.58. For 
certain values the probability of detecting photons can 
be larger for balanced functions than for constant ones, 
so we would infer the function wrongly more often than 
correctly. Of course, we would then have to bet in the 
opposite way: balanced if we detect a photon, constant 
if we do not detect anything. For a detector covering the 
whole plane we recover the classical fifty-fifty guess since, 
again, we do not get any information about /. 



We note finally that in this setting we have an asym- 
metry between the betting confidence on constant and 
balanced functions: For small detector sizes, a detection 
implies a constant function with high probability, while 
no detection implies a balanced one just with moderate 
probability. Indeed, from Bayes' formula we can com- 
pute: 

P(f = constant | detection) = — , (7) 

Pc +Pb 

1 — ^iPb 

P(f = balanced I no-detection) = ; r-. (8) 

1 ' ' 2- V { Pc + Pb ) y ' 

For instance, for a very small detector we have p b <C 
p c < 1 so P(f = constant | detection) « 1 while P(f = 
balanced | no-detection) « 1/2. This is due to our choice 
of the detector position. Choosing a spot on the Fourier 
plane corresponding to the state |— ) would invert this 
asymmetry. If a CCD was used instead we would have, 
on average, the same confidence for betting on both 
types of functions, but there would still be inconclusive 
events. If the signal photon were detected by the CCD 
on a position corresponding to the state preparation, say, 
^(|0) +i |1)), we could not infer the type of function we 
had. For instance, a detection at position O-llrnm, i.e., 
where the two interference patterns cross in Fig. [2j cor- 
responds to such a situation. 
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FIG. 4. Probability of success when the signal detector size 
is changed (width of the single slit in front of the detector is 
varied ). The photon detectors are assumed to be fixed at 
the positions Xi — x a — 0, and it is assumed the detection of 
the idler photon. Here the detectors are considered perfect, 
and we consider the average probability of both detection and 
no-detection events. It is assumed that we always bet in a 
constant (balanced) function when there is (no) detection. 
The optimal size for the detector is 260/im, with a success 
probability of 0.58. 
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VI. CONCLUSIONS 

In this work we implement the minimal Deutsch algo- 
rithm version with one qubit, using spontaneous para- 
metric down conversion and idler photodctection as a 
heralded source of one photon and a spatial light modu- 
lator as the key part of the quantum oracle. A double- 
slit is used to encode a qubit in the photonic transver- 
sal spatial modes and the state is manipulated using the 
SLM. The experimental setup is able to implement all 
possible one bit constant and balanced functions easily. 
Furthermore, we discuss the improvement of this specific 
quantum experimental implementation, when compared 
to the analogous classical algorithm. This is the first 
quantum algorithm implemented using a SLM, opens the 



possibility to implement more complex algorithms. For 
example, by changing the double-slit to an eight-slit in- 
terferometer [l5j], one can also implement the analogous 
minimal version of Deutsch- Jozsa algorithm for three-bit 
functions / : {0, l} 3 -> {0, 1}. 
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